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Abstract. Star formation depends on the available gaseous "fuel" as well as galactic environ- 
ment, with higher specific star formation rates where gas is predominantly molecular and where 
stellar (and dark matter) densities are higher. The partition of gas into different thermal com- 
ponents must itself depend on the star formation rate, since a steady state distribution requires 
a balance between heating (largely from stellar UV for the atomic component) and cooling. 
In this presentation, I discuss a simple thermal and dynamical equilibrium model for the star 
formation rate in disk galaxies, where the basic inputs are the total surface density of gas and 
the volume density of stars and dark matter, averaged over ~ kpc scales. Galactic environment 
is important because the vertical gravity of the stars and dark matter compress gas toward the 
midplane, helping to establish the pressure, and hence the cooling rate. In equilibrium, the star 
formation rate must evolve until the gas heating rate is high enough to balance this cooling rate 
and maintain the pressure imposed by the local gravitational field. In addition to discussing the 
formulation of this equilibrium model, I review the current status of numerical simulations of 
multiphase disks, focusing on measurements of quantities that characterize the mean properties 
of the diffuse ISM. Based on simulations, turbulence levels in the diffuse ISM appear relatively 
insensitive to local disk conditions and energetic driving rates, consistent with observations. It 
remains to be determined, both from observations and simulations, how mass exchange processes 
control the ratio of cold-to-warm gas in the atomic ISM. 



1. Introduction 

Disk galaxies are gas-rich systems, with a multi-phase, highly structured interstellar 
medium (ISM). Within the ISM, star formation takes place in giant molecular clouds 
(GMCs), sometimes concentrated in spiral arms. The rate and charact er of star formation 



are in fluenced by physical processes from sub-pc to multi-kpc scales (jMcKee fc Ostriker 
2007). In spite of the complexity of the ISM and star formation at small scales, there are 
nevertheless clear correlations between the large-scale rate at which stars are born, and 
the properties of the ISM and (intra-) galactic environ ment on large (~ k pc) scales. 
As discussed by Frank Bigiel at this meeting (see also Bigiel et al.ll2008l and references 



therein), in regions of galaxies where the gaseous surface density E £ 100 M Q pc~ 2 , 
the star formation rate closely follows the surface density of molecular gas. This can be 
understood in terms of having an essentially constant star formation timescale, t^F ~ 
2 x 10 9 yr, within molecular gas (which is observed to be in organized in gravitationally 
bound clouds with properties that are similar in different galaxies). As a consequence, 
Ssfr oc E in regions where the molecular gas dominates the atomic gas. For regions where 
atomic gas dominates (primarily in the outer parts of galaxies) , Esfr instead varies as a 
steeper power of E. In addition to this superlinear behavior, there is considerable scatter 
in the relation between Esfr vs. E at low surface density, suggesting that one or more 
other parameters, in addition to E, controls the star formation rate. 
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Indeed, recent examination of the correlation of Esfr with "non-interstellar" galactic 
environmental properties has revealed interesting dependences, indicating that in the 
outer parts of galaxies, both the specific star formation rate and the ratio o f molecular 
to-at omic gas incr ease roughly linearly with t he stellar surface density E s (jLerov et al 



20081) . Previously, iBlitz fc Rosolowskvl (j2006l ) found an approximately linear increase 
of the molecular content with the estimated dynamic pressure of the ISM, and this is 
evident in the sample analyzed bv iLerov et al.l (21)08) as well. The physical reason for 
the relationship between molecular content (and star formation) and pressure has not, 
however, been clear from these empirical studies. 

Observations of star formation pose a number of challenges: Why is there an increase in 
the slope of Esfr oc E 1+p in going from molecular- to atomic-dominated regions? What 
is the physical reason for the empirical relation between ISM pressure and star formation; 
more generally, how do galactic parameters such as E s , the velocity dispersions of stars 
and of gas, and spiral structure affect Esfr? Is it possible to explain the observed behavior 
of Esfr using simplified theoretical models, and what is required in numerical simulations 
in order to reproduce observed star formation relationships? Recent theoretical work has 
taken on these challenges with increasing success; a key to these advances has been a 
more sophisticated t r eatmen t of both the ISM and the galactic environment. For example, 
Kovama fc Ostriker ( 2009al ) found, using numerical simulations of the ISM and a cooling 



function allowing multiple phases, star formation rates and proportions between self- 
gravitating and di ffuse gas similar to the observations of Blitz fc Rosolowskv ( 20061 ) and 
Lerov et al. I (120081) provided that turbulent driving is included; for non-turbulent models, 



the proportion of self-gravitating gas was found to be much too high. 



2. A thermal/dynamical equilibrium model for Ssfr 

Mot ivated by recent observations as well as simulations and earlier theory. lOstriker et al 



( 2010f) (hereafter OML) have developed a simple model for star formation regulation in 
multiphase, turbulent ISM disks. In essence, the OML model combines three basic prin- 
ciples: (1) the diffuse (atomic) component of the ISM is in approximate thermal equi- 
librium, with a density (and pressure) proportional to the heating rate; (2) the diffuse 
component of the ISM is in approximate dynamical equilibrium, with the pressure at 
any height above the galactic midplane given by the weight of the overlying gas; (3) 
UV from young stars provides most of the heating for the atomic component of the 
ISM, with star formation taking place only within the gravitationally-bound component 
of the ISM. These principles have been individually es tablished and exte nsively studied 
(over several decades) in the astrophysical literature. Field et al] ( 19691) combined (1) 
and (2) to conclude that the diffuse atomic gas in the local Milky Way must consist of a 
two-phase cloud/intercloud medium. In this and subsequent treatments of thermal and 
dynamical equilibrium, the heating rate has generally been treated as an independent 
(empirical) parameter. But, by including (3) together with (1) and (2), OML obtained 
a closed system representing a local patch of a disk galaxy. For this closed system, the 
partition of gas into phases and the star formation rate are obtained self-consistently. 

In the OML model, the (simplified) ISM is treated as having two components, one 
consisting of diffuse gas (including both high-density cold atomic cloudlets and a low- 
density warm atomic intercloud medium), and the other consisting of gravitationally- 
bound clouds (GBCs). Although hot gas is also present in the ISM, it is a tiny fraction 
of the mass, and fills ,$ 20% of the volume ( Heilesll200ll ) (OML describe how to correct 
for this effect). For galaxies with normal metallicity, the GBCs would represent giant 
molecular clouds, including their atomic shielding layers. Averaged over ~ kpc scales 
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(which may contain many or few individual GBCs), the total surface density of the GBC 
component is Sgbc, and the total surface density of the diffuse component is Ediff- 

The diffuse component is assumed to be in vertical dynamical equilibrium (as has 
been v erified by numerical simulations; e.g. Piontek fc Ostrikerll2007 ; Kovama fc Ostriker 



l2009bl ). with the vertical gravity (from the diffuse gas, the GBC component, the stellar 
disk, and the dark matter halo) balanced by the difference between midplanc and external 
values of thermal pressure P t h, turbulent pressure pv%, and magnetic stresses (87r) _1 (B 2 — 
2B 2 ). Because cooling times are short compared to other timescales, the diffuse gas is 
assumed to be in thermal equilibrium, with the additional provision that both warm and 
cold phases are present. This allows a range of pressures between P m i n ,coid and P m ax,warm; 
for the model of OML, it is assumed that the pressure is equal to the geometric mean of 
these limits, Ptwo-phasc = (Pmin,coidPmax,warm) 1/2 - For atomic gas, heating is generall y 
dominated by the UV and cooling by collisionally-excited lines ( Wolfire et al.l 120031 ) . 



which yields Ptwo-phasc oc Juv- (Note that other heating - e.g. cosmic rays and shocks 
- can be more important for very dense, shielded cores and very hot gas, respectively.) 
Finally, the OML model assumes that the rate of star formation is proportional to the 
total surface density Sqbc of gas in the GBC component, Esfr = SqbcAsf = (S — 

£diff)AsF- 

Vertical dynamical equilibrium within the diffuse layer is expressed as Ptot = oPth = 
Sdiff(gz)/2, where the mean vertical gravity is 

(g z ) « ^G(S diff + 2S GBC ) + 2(2G (0sd ) 1/ V z ; (2.1) 

p s d is the midplane density of stars plus dark matter, o~ z is the total vertical velocity 
dispersion of the diffuse gas, and the total pressure is larger than the thermal pressure 
by a factor a (see below). The GBC component contributes more strongly (per unit 
mass) to the gravity because its scale height is smaller than that of the diffuse gas. 

If n 2 A(T) is the cooling rate per unit volume and nT is the heating rate per unit 
volume, then the two-phase pressure is given by 

Ptwo-phasc i rp rp 

T = V^min, cold-' min,cold^'max, warm-' max, warm ) 

p (Jmin, cold Jmax, warm) ^ 

[A (Jmin, cold) A(Tn± ax, warm)] 

where we have used the equilibrium condition T = nA for both phases. Cooling of the 
cold atomic medium is dominated by metals (in particular, C II) so that A oc Z gas , 
while heating is dominated by the photoelectric effect with T oc ^dust^uv; since r m in.coid 
and Tmax^arm are relatively independent of the heating rate ( Wolfire et al. 19951 ). this 



yields Pt wo -phasc oc Juv if Zd us t/Z gas = const. The terms Z gas and Zdust represent the 
ratios of metals and dust to hydrogen, respectively. The mean UV intensity is affected 
by radiative transfer through the diffuse gas, but for modest optical depth in the diffuse 
gas the relation Juv Ssfr is expected to hold. In addition, a larger fraction of the UV 
escapes from GBCs if Zd is very sub-Solar, which increases Juv for a given Esfr (this 
effect is quite uncertain, but might increase Juv by a factor ~ 2). In thermal equilibrium 
with P tri ~ Ptwo-phasc, the midplane pressure is therefore expected to vary roughly as 
Pth oc Ssfr, with a somewhat larger coefficient for very low-metallicity regions. 

Combining the thermal equilibrium relation P t h = Pth,o^SFR/SsFR,o (normalized us- 
ing the Solar neighborhood thermal pressure Pth.o and star formation rate Ssfr,o) with 
the dynamical equilibrium relation P t h = Sdiff (9z)/(^ a ) and the star formation relation 
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E SF r = EgbcAsf, we obtain 

^£2£ = M a ^G(E diff + 2E G bc) + 2(2G Psd ) 1 / 2 a 2 . (2.3) 
^diff g* 

Here, = 2aPth,o/ (^sfr.o^Sf); for fiducial parameters, this acceleration is <?* ~ pc Myr~ 2 . 

It is interesting to compare outer and inner disks. In outer disks (similar to the Solar 
neighborhood and beyond, in galaxies like the Milky- Way), diffuse gas dominates the 
total so that Eqbc *C E d ig ~ E; in addition, the term depending on /? sd dominates the 
gravity g z . In this regime, the relation Esfr oc Eqbc ck ^yffhd is therefore expected to 
hold. Physically, this regime may be thought of as the result of star formation increasing 
until the heating it provides is sufficient to balance cooling at the (dynamically-imposed) 
midplane pressure. If there is too little gas in the GBC component, the star formation rate 
would be extremely low, and the UV field would be very weak. A very low heating rate 
could not maintain a warm medium at the pressure imposed by the local gravitational 
field, so that (a portion of the) warm gas would condense out and become cold clouds. 
These cold clouds would collect to create more GBCs, which would then initiate star 
formation, raising the local UV radiation field until heating balances cooling. Given the 
low gravity and pressure of outer disks, cooling rates are moderate, and relatively low 
levels of star formation are needed to produce enough UV that heating balances cooling. 

For outer disks where the stars and dark matter dominate gravity, the vertical oscilla- 
tion time is t osc — 7r 1 / 2 /(G j o s d) 1 ^ 2 ; a dense cloud settles to the midplane in ~ t osc /4. In 
this regime, the conversion time from gas to stars, t con = E/Esfr, is given by 

* -t azPtK0 (2 4) 

'con — <-osc ,„ M/2/ 2 \v ' V^-^V 

O") 1 2 Kh) S SFR,0 

where (w t 2 h ) = ] w (? w w c^M d i ff)Warro /M d i ff)tota i is the mean thermal dispersion in the 
diffuse medium (here c w ~ 8 km s" 1 is the thermal speed in the warm ISM). Using 
-Pth.o ~ (^h^gas.o/ " 2 an d defining a star formation energy conversion efficiency e rad = 
47rJ rad>0 /(c 2 I]sFR,o) for -Prad.O — 47rJ radi0 /(3c), 

icon = tosc 7777, vT72 D £ rad- (2-5) 

3(27r) i /^ - z P iadtQ 

That is, the gas conversion time (or depletion time) is set by the time for gas to settle 
to the midplane, scaled by factors for the ratio of gas-to-radiation pressure in the Solar 
neighborhood, the mass-to-energy conversion efficiency, and c/a z . 

In inner disks, unlike outer disks, we have S d ig <Si Sgbc ~ so that Esfr oc S. 
In inner disks, it is straightforward to show that there is an upper limit on the diffuse 
gas surface density S d iff. Physically, the reason for this limit is that the diffuse-gas 
cooling rate per particle increases with higher density and pressure in the inner parts of 
disks at least as nA oc S d ;ffSGBC (since nA oc S d ;ff/iJ oc Ediff^/c 2 oc E d iffEGBc[l + 
ffsd/ffGBc]/c 2 )i whereas the heating rate per particle varies as T oc Ssfr oc Eqbc- Thus, 
cooling will exceed heating (causing mass to drop out of the diffuse component) unless 
Ediff is sufficiently low. Enhanced cooling and mass "dropout" is likely responsible at 
least in part for the "saturation" of HI sur face densities at ^ 10 M Q pc~ 2 that has 
been observed in the inner parts of galaxies ( Bigiel et al. 20081) . 



Based on the relations described above, the star formation law is expected to steepen 
from Esfr oc E in inner disks to Esfr oc E^/p sd in outer disks. A reduction of the specific 
star formation rate Esfr /E is indeed observed in galaxies starting at E ,$ 10 M Q pc -2 
(|Bigiel et all 120081 : iLerov et al.ll2008l ). For some galaxies, a further power-law relation 
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Figure 1. Comp arison between annular averages of the data (squares) for NGC 5055 
ijLerov et alj200a i. and the thermal/dynamical equilibrium model (triangles) developed in OML. 
Both the star formation rates as a function of radius in the galaxy (panel d), and star formation 
rates as a function of gas and stellar density (panels a and b) agree with the model predictions. 

p s d oc £ 2p may hold such that Ssfr oc in outer disks, but this need not be the case 
in general - that is, integrated "Schmidt" -type relations need not hold. 

In OML, the full solution for Ssfr is obtained as a function of E, p S( j, and the param- 
eters a = Ptot/-fth and f w = ( v th)/ c wi under the assumptions of thermal and dynamical 
equilibrium described above. It is also shown that the theoretical solution for Ssfr agrees 
well overall with a sample of disk galaxies analyzed in lLerov et al. (2008), with especially 



close correspondence for the large flocculent galaxies NGC 7331 and NGC 5055. Figure 
1 shows an example of the comparison between the model and data, for NGC 5055. 

Given the promising comparisons between the analytic theory and observations, it will 
be quite interesting to develop numerical simulations that fully test the assumptions and 
results of the thermal/dynamical equilibrium model. Encouragingly, the poster presented 
by C.-G. Kim at this meeting shows that initial numerical tests support the assumptions 
of thermal and dynamical equilibrium adopted in the analysis of OML. As discussed 
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above, the OML theory contains parameters that must be set from either observations 
or detailed simulations. In the remainder of this contribution, we review what is known 
in this regard based on previous numerical work, and what measurements will be needed 
from future modeling efforts. 



3. Numerical evaluation of parameters 

From equations (12.11) and (|2.3p , the star formation rate in outer-disk regions is expected 
to vary as Ssfr oc ^y / 2Gp sd cr z /a, where a = a 2 /v 2 h and a 2 = v% h + v% mh + (l/2)A(v\ - 



ith 



'th' ^turb' 



and v A the (mass- weighted) mean thermal, turbulent, and Alfvcn 
speeds in the diffuse gas (we now omit angle brackets denoting averaging) . The coefficient 
<T z /a can also be written as v\J a z — c^jwj o z . Thus, the star formation rate is expected 
to depend on the total velocity dispersion a z (or the ratio u z /c Wl where c w is fixed by 
atomic physics), and on the fraction of diffuse gas in the warm phase ss f w 

The ratios a z /c w and f w Mdiff,warm/Miiff,totai depend on the details of gas dynamics 
in the diffuse ISM. Important effects include warm and cold phase exchange via thermal 
instability; turbulence (with the associated shock heating and adiabatic temperature 
changes, as well as turbulent mixing); conversion of diffuse gas to GBCs via midplane 
settling, self-gravity, and turbulence-induced cloudlet collisions; return of GBC gas to the 
diffuse phase by photodissociation and by "mechanical" destruction processes (including 
expanding HII regions, winds, SNe, and radiation pressure). Turbulence in the diffuse gas 
can be driven by stellar energetic inputs as well as spiral shocks, the magnetorotational 
instability, large-scale gravitational instabilities in the disk, and cosmic infall. 

Numerical studies to understand the various effects involved are very much a work in 
progress, but some consensus is already beginning to emerge on a number of points: 

• For a medium with a bistable cooling curve, the midplane thermal pressure tends to 
evolve, by exchange of mass between cold and warm components of th e diffuse phase, such 
that t he me an value is comparable to, or slightly below, Ptwo-phase (jPiontek fc Ostriker 
200.4 l2007h . Since out-of equilibrium effects depend on the heating time from shocks 
compared to the cooling time, the mean value of the thermal pressure, as well as the 
breadth of the pressure distri bution, must in general be affected by the scale and the 



amplitude of turbulence (see iGazol et al.l 12003. 120091: lAudit fc HennebelleJ 120051 lioiot 



Hennebelle fc Auditl[2007t I Joung fc Mac Lowll2006l : I Joung et al.ll2009f ). Realistic numeri 



cal evaluations of the mean thermal pressure (for a given radiative heating rate) therefore 
will require numerical simulations in which the vertical box size is comparable to the true 
scale height of the diffuse ISM, and in which the turbulent amplitude is ~ 5 — 10 km s -1 . 

• Magnetic fields in differentially-rotating multiphase disks are amplified by the mag- 
netorotational instability until the magnetic pressure becomes comparable to the ther - 
mal gas pressure, with B 2 < B 2 (jPiontek fc Ostriker|[200l l2007t I Wang fc Abe3l2009[). 



Supernova-driven turbulenc e also contributes to amplifying the magnetic field (Ide Avillez fc Breitschwerdt 
20051 : iMac Low et alJliooil ). 

• The energy input from supernovae yield ISM velocity dispersions ~ 5 — 10 km s _1 for 



models with a wide range of supernova driving ra t es and disk properties (e . g. lKimll2004 



de Avillez fc Breitschwerdtl 120051 : iDib et al.|[200l IShettv fc Ostriker|[20oj lAgertz et al 
20091 : 1 Joung et al.ll2009h ~ These values are comparable to those observed in the HI gas 
Simulations have also shown that the turbulent amplitudes decrease at smaller scales 
and for higher densities. With this range of turbulent velocity dispersions, the turbulent 
pressure in simulations of the diffuse ISM is comparable to the thermal pressure. 

• The interaction between self-gravity and ro tational shear also drives turbulence at 
significant levels ( ^ 10 km s _1 ) in galactic disks ( Kim fc Ostrikerl2001 . 2007 : Wada et al 
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2002tlShettv fc Ostrikerll2008t iTasker fc Ta"n1l2009l; lAgertz et al]|2009t lAumer et allbOld; 
Bournaud et al. 2010h . However, the turbulent power is much larger at the large (~ kpc) 



scales that dominate the swing amplifier than at scales below the disk scale height, and 
in-plane velocities (which do not contribute to vertical support of the disk) are much 
larger than vertical velocities. Thus, turbulence driven by instabilities on large scales is 
likely of limited importance in regulating the effective midplane pressure (for a given lo- 
cal gas surface density E), and hence the star formation rate. (Gravitational instabilities 
would, however, enhance S and thus Esfr locally.) Fl apping associated with non-stead y 
spiral shocks also drives turbulence in the diffuse ISM (|Kim et al.ll2006 ; Kim et al~ll2010[ ) , 
but again, vertical motions are small compared to horizontal motions. 

Although numerical results have shown that the total turbulent velocity dispersion o~ z 
is relatively insensitive to the disk properties and the supernova driving rate (consistent 
with observations), it is much less certain how the warm fraction, or v^ h = j w c 2 w ~ 
c^Mdiff jWa rm/-Wdiff, total, depends on disk conditions and/or the star formation rate. As- 
sessing this dependence, including a full exploration of parameter space, is an important 
task for future numerical studies. The fraction of diffuse atomic gas in different phases 
is not well known empirically, either, although observations of C II with Herschel po- 
tentially afford a means to separate cold and warm components of the atomic medium 
(which both contribute to 21 cm emission). 

Finally, it remains important to understand more fully how spiral structure develops, 
and in particular, whether it is possible to characterize in a simple way the fraction 
of gas in a given annulus that is found in "arm" vs. "interarm" conditions, and what 
the compression factor is for the gas surface density. Numerical simulations have begun 
to marry spiral structure with an increasingly realistic treatment of the ISM (including 
multiple phases, turbulence, and magnetic fields); much more, however, remains to be 
done on this front. It also remains to be determined how well models like that of OML 
apply locally for galaxies with strong spiral structure. More generally, it is important 
to assess which equilibria (thermal, dynamical, star formation) still apply locally even 
in galaxies with large-scale transient structure in the ISM (due to spiral arms, tidal 
interactions, mergers, cosmic inflows, etc.). 



4. Conclusion 

Gas is the raw material for star formation, but the detailed state of the ISM, which 
depends in turn on the internal galactic environment, determines the rate at which this 
material is processed to create new stars. Recent observations have begun to explore the 
correlation between gas content and star formation at increasingly high spatial resolu- 
tion, revealing changes in star formation "laws" between inner and outer disks; other 
environmental dependences of star formation have also been explored, including intrigu- 
ing correlations between molecular and stellar content of galactic disks. 

Although the simplest recipes for star formation (such as a rate that depends inversely 
on the free-fall time at the mean ISM density) have difficulty matching the data, models 
that account for feedback and the multiphase character of the ISM are more successful. In 
particular, recent work suggests that the empirical correlation between molecular content 
and estimated midplane pressure can be understood as reflecting a state of simultaneous 
thermal and dynamical equilibrium in the diffuse ISM. The thermal/dynamical equilib- 
rium model of OML develops the idea that UV from OB stars provides a feedback loop 
that regulates the star formation rate: the proportions of diffuse and self-gravitating gas 
adjust themselves so that the heating rate (proportional to the mass of self-gravitating 
gas) matches the cooling rate (proportional to the mass of diffuse gas and to the vertical 
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gravitational field). The model formulated in OML is promising in terms of explain- 
ing star-forming behavior in observed systems. With numerical simulations, it will be 
possible to appraise - and potentially revise - the simplifying assumptions and param- 
eterizations adopted by this equilibrium model. Time-dependent simulations will also 
lead to a much clearer understanding of how GBCs form and disperse, and how their 
properties and the formation/destruction timescales relate to galactic environment. This 
will aid in defining limits for applying equilibrium relations, while also pointing the way 
towards non-equilibrium theories of star formation. 
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